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ABSTRACT 

For each photon wave packet of extragalactic light, the dispersion by line-of- 
sight intergalactic plasma causes an increase in the envelope width and a chirp 
(drift) in the carrier frequency. It is shown that for continuous emission of many 
temporally overlapping wave packets with random epoch phases, such as quasars 
in the radio band, this in turn leads to quasi-periodic variations in the inten- 
sity of the arriving light on timescales between the coherence time (defined as 
the reciprocal of the bandwidth of frequency selection, taken here as of order 
0.01 GHz for radio observations) and the stretched envelope, with most of the 
fluctuation power on the latter scale which is typically in the millisecond range 
for intergalactic dispersion. Thus, by monitoring quasar light curves on such 
short scales, it should be possible to determine the line-of-sight plasma column 
along the many directions and distances to the various quasars, affording one a 
3-dimensional picture of the ionized baryons in the near universe. 



Introduction 



After the reionization epoch, the diffuse baryons of the universe has a plasma component 
that for z < 1 redshifts is predominantly within the 'warm' temperature range 10 5 ~ 6 K, and 
accou nting for 40 - 50 % of the baryons in the near universe (ICen and Ostriker I (119991 ) ; iDave 



(120011 )). This form of matter fills the entire intergalactic medium (IGM) filament arily, with 
concentrations near clusters and groups of galaxies. It still largely evades detection because 
the emission is in the EUV and soft X-rays, a wavelength passband that suffers from Galactic 
absorpti on as well a s inter stellar and time variable heliospheric foreground contamination, 



(see e.g. 



Takei et al I (120081) on the last; and fo r up to date inform ation on the search o f these 



'missing' baryons see iBonamente et al I (120121 ) and the review of iDurret et al I ( 120081 )). For 
this reason other ingenious and 'orthogonal' techniques, including and especially the recent 
search b y the Plan c k team of 'excess' Sunyaev-Zel'dovich (SZ) signals at the outskirts of 
clusters (Ude et al I (l2012al lbh may be very important, as both papers reported positive 
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results. In par ticular, t 



Coma cluster (lAde et ai- 



r e 'exce ss' Compton y-parameter at the 2-5 Mpc radii of the 
(|2012a|)l is consistent wit h the predictions based upon soft X-ray 



observations of Coma ( Lieu and Bonamente ( 20091 )). as both SZ and X-ray data yielded the 
same 5 Mpc limiting radius for the extent of the warm baryons as well. 

Nonetheless, the warm baryons are not necessarily confined to clusters' vicinity alone, 
but can also structurally fill the rest of the IGM so that even the SZ technique is unable 
to deliver a complete search. This is because the technique suits regions of high baryonic 
density and temperature, viz. baryons inside clusters can exert a large SZ 'pressure' along 
the line-of-sight of interest, whereas outside the clusters the plasma becomes cooler and more 
tenuous, so that this 'pressure' can easily fall below SZ detection thresholds. The search for 
genuinely intergalactic baryons must therefore await new approaches. 

Here we intend to discuss one of them. Although the primary difficulty with the po- 
tentially very powerful method of utilizing the light from quasars (the only point sources 
bright enough to be visible across cosmological scales) as probes of the i ntervening un i- 



verse is the lack of fast pulsar like i ntrinsic variations in their intensity (ILazio I (120081 ); 



Dennett-Thorpe and de Bruyn I (120021 ). the latter showed that rapid (hourly or shorter) 
variabilities in the quasar light are not intrinsic to the source), it will be argued in this 
Letter that the plasma dispersion effect of photon envelope broadening and carrier wave 
chirping are imprinted upon even steady and continuous light signals as they pass through 
the vast spans of the IGM. The imprint can be uncovered by conventional techniques if 
the observations are done at radio frequencies. Previous effort in this vein were mainly 
done with n on-astrophysical applicatio ns in mind. They include measurements^ of coher- 



ence length (IHitzenberger et al I (jl999|)), and calculation of the broadening of wide pulses 
( Saleh &: Irshid 1 ( 1982 )) and pulse distortion to 3rd order Taylor series correction ( Marcuse. D. 
(Il980l )) in a dispersive medium. To the best of the authors' knowledge there has not been 
any treatment of the time resolved mutual coherence function of continuous light and its 
associated noise characteristics in the same. 



2. The dispersion of quasar light 

We begin by revisiting the question of how, in an astrophysical context, a single wave 
packet of light behaves as it passes through a dispersive medium. Since the cosmological 
journey of light from an unresolvable point source to a small telescope on earth is principally 
a one dimensional problem k = (k,0,0) one can under this scenario ignore the dynamics of 



x For the theory see http://light.ece.illinois.edu/ECE460/PDF/LCI.pdf. 
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the wave packet in the y and z directions, i. e. the spreading of the packet is only appreciable 
along x (scattering by plasma clumps can cause broadenin g, but for a mean p lasma column 
density equivalent to 1 Gpc the effect is at the 10~ 6 s level (IBhat et al I (120041 )). which as we 
shall see in section 3 is ~ 10 4 times below that of plasma dispersion). The amplitude of a 
1-D gaussian wave packet emitted at t = t e (more precisely the packet's center is x = x e at 
t = t e ) at a point (t,x) 'downstream', may be written as 



ip{t,x) 



.Ma) 



2(Afc)2 e i[k(x-X-)-U)(t-t e )] flfc 



(1) 



where the spectral filter is oc \f(k)\ 2 with 

/(*) 



g 2(Afe) 2 



(2) 



and is peaked at k$ of width Ak. The filter / usually depicts either line emission at the 
source or passband selection by the observer, the two are equivalent because the plasma 
medium does not distort the frequency spectrum of the radiation. 

An approximate form of (pQ) that reveals its salient features is afforded by Taylor ex- 
panding u = oj(k) around ko to 2nd order (3rd and higher order terms can be ignored in the 
case of plasma dispersion provided Ak/k <C 1), (pQ) then becomes 



if}(t - t e ,x) 



A 



2w(Akf 



l + iuj{(Ak) 2 (t-t e 



exp 



(AA;) 2 [x - x t 



exp [ik (x — x e ) — iuo(t — t e )\ 
v g (t-t e )} 2 [l-t^(t-t e )(Ak) 2 ] 



where y s = {<ko/ dk) k=k n and u'q 



m e.g. 



2 l + uol\AkY{t-t e ) 2 

(d 2 uj/dk 2 )k=k - Details of the derivation of 



(3) 



are given 



Bohm I (jl95fl ). section 3.5. 



We turn to look at the nature of continuous light from a cosmological source that is 
steady on short time scales. The amplitude of which at some instance t and position r, 
with t as a variable and r held fixed, may be expressed as the sum of many overlapping 
pulses of amplitudes ipj(t), each depicting a single photon consisting initially of a harmonic 
(i.e. unchirped) carrier wave with random phase and enveloped by the coherence time r c , 
viz. i.e. 

V>(*) = X>*(*)> andr c = — . (4) 



In our present labeling scheme the photon pulses are chronologically ordered in terms of 
their envelope peaks. The actual number n of such pulses that significantly influence any 
given instance t is finite, however, due to the finite width r c . This is the situation before 
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dispersion. In fact, n is the 'occupation number', or the number of arriving photons per 
bandwidth per coherence time that becomes an invariant for a given source and telescope. 

The amplitude function of the arriving light, however, may be quite different. It is still 
given by of course, but the constituent wave function ipj may now be modified, viz. 



A / 27T 

= — J 7^ + i£) Au ex P[ _a (* _ ^) 2 ] ex P[ ife (* ~ ^i) 2 ~ iu} o(t - tj) + <f>j], (5) 

where 

a = ^>n~T~c2\' and b = ^ Wlth £ = — 2 ' ( 6 ) 

2(1 + 4) 

and 0j a random phase§. Note that £ changes very slowly with time and may be treated 
as a constant. In an expanding universe, it is determined by the comoving distance to the 
source, and the carrier frequency and average IGM plasma density (with both evaluated at 
the present epoch), as we shall find out. 

Let us first examine the microscopic variability of the arriving intensity / = ^(t)] 2 
which may be broken down into two parts, 1 = 1 + 1%, where 

n'+l n'+l 

i = E i^wr> and h = E (?) 

With the help of (JSJ) , one may in turn write I as 

2 p_ 2 e _ 2a(t _ tj)2 = |A|2 (Au;) 2 = \A\ 2 V^(Aw) 2 , (8) 

1 1 C 2 v /rTe 2V cV2(i + e) 1 1 V2c 2V ; 



where the rightmost expression applies to the limit of large occupation number n>l (hence 
n' ^> 1 necessarily) when 7 becomes a constant. 

Next we work on 1%, which appears as 

n'+l 



h = h(t) = \A\ 2 -^===(Au;f e^-^e-^-^ 2 cos[26(t fc - + <p jk ], (9) 



1 v ' ^ i/fc 



2 Our treatment here does not apply to the scenario of phase coheren c e amo ng many temporally contiguous 
photons, i.e. the phenomenon of photon bunching (e.g. iRickett et al I (|1975f )) which is unlikely to apply to 
intrinsically slow emitters like quasars. 
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where cfjk = —<Pkj changes randomly from one distinct pair of j, k to the next. The mean of 
I\ obviously vanishes. The variance is given by 

9-7T 2 n ' + 1 

a] = |A| 2 -^— (Ac) 4 £ e- 2 ^) 2 e--^) 2 (cos 2 [26(4 - *,■)* + <p jh )). (10) 

Again, in the large n limit the double sum may be approximated by an integral, and (cos 2 ) = 
1/2. The result is 

Gl = \ A \ 2 (Ac) 2 = \A\ 2 ^^(Au) 2 = -. (11) 

2cV2(l+e) 2V2c 2K ' 2 v ; 

Thus, both / and 07 are independent of dispersion - they are only functions of the photon 
occupation number n and bandwidth Ac (or coherence time r c via (jib)). 

The interesting question is the timescale over which the intensity I undergoes the random 
variation about the mean / with the variance of (11 II) . The answer comes from the cos[2o(£fc — 
tj)t + ipj k ] factor of Under the scenario of no dispersion, i.e. a = (Ac) 2 /2 and 6 = 0, 
the factor only depends on the set of phase angles {fjk}, which changes to a completely 
different set when the time t translates by ~ one envelope width 1/Ac, or r c . Thus, in 



accordance with known facts (IMathieu. J.P. I (Il975[ )). there is a quasi-periodic variation in 



the intensity on the scale of one coherence time, and the Fourier Transform of / would show 
a spike at the frequency ~ Au. 

After passage through dispersion, however, the cosine factor will still vary significantly 
at the level of ( ITTj) when {(fijk} is replaced by a completely different set, except this now 
happens on the timescale of one envelope width a -1 / 2 ~ t c a/2(1 + £ 2 ) 3> r c ; moreover there 
will also be milder fluctuations over shorter timescales due to the 2b(tk — tj)t part of the 
cosine argument. To elaborate, in Yljk ^ ne Q uan tity — tj ranges from ~ r c (the assumed 
timing accuracy of intensity measurements) to ~ a~ 1//2 . In the former end of the range 
26(tfe — tj)5t — > 1 when 5t ~ a -1 / 2 , while in the latter end 5t ~ r c . Since each equal interval 
of tk — tj has about the same number of photons, this means fluctuations at the level of a 
fraction of ( TTTT) exist on all scales between the coherence time r c and the envelope width 
a -1 / 2 , with the Fourier power increasing towards the scale of a -1 / 2 where it reaches the full 
level of (TTTT) because of the behavior of {<Pjk}- This is the observable imprint of dispersion 
upon the passing light: the microscopic details of the intensity variation are very different 
from the scenario of a non-dispersive medium where only the characteristic timescale of r c 
exists. 



We end this section by presenting a self-consistency check of the formalism of ([7]) and 
05]). They can also be used to demonstrate a known fact: the independence of the coherence 
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time of continuous light, denned as the maximum time delay for observable interference 
fringe contrasts, on dispersion (so long as the medium does not distort the radiation energy 
spectrum). The proof involves calculating the autocorrelation function (ACF) at the delay 
r, or 2 J dt Re[if}(t)ijj*(t + r)], from the two equations in much the same manner as above, 
to obtain the result 

2 J dt Re[ip(t)ip*{t + r)] = Icosuot e~ r2/{ ^\ (12) 
which shows that the periodic fringe pattern cos coqt is damped exponentially away when the 



delay r exceeds r c (the van Cittert-Zernike theorem, see e.g. section 10.4.2 of iBorn fc Wolf 



( 1970 )). irrespective of dispersion^ as £ does not appear anywhere here. Moreover, because 



in (TT2"j) the pattern as given by the ratio of the ACF to / is not a function of n, each photon 
interferes only with itself. 



3. Measurement of IGM plasma column density 

For dispersion in a cold plasma the dimensionless quantity in (j6b) is given by 



u \ 1 / n e 



(13) 



6 x 10 9 rad/s J V lO" 7 cm.- 3 J \1 Gpc 

The central received frequency of Vq = 1 GHz would correspond to radio observations. In (fT3"l) 
we ignored a small relativistic correction due to the fact that the IGM plasma at redshifts 
z < 1 is warm. Note also that in f fl3|) we ignored the expansion of the universe, which does 
not introduce significant errors for sources with £ < 1 Gpc. If account is taken of the in situ 
expansion, the dimensionless parameter £ will become 

^f^'^i (U) 



^0 / w J te a(t') 

where u p and uo are respectively the IGM plasma frequency and the radiation frequency as 
evaluated at the time of observation t, t e is the time of emission, and a(t) is the expansion 
parameter. From (jT4"|) it becomes apparent that, although the value of £ for sources at large 
distances require redshift correction, this is made simply by interpreting I in ffl3|) as the 
comoving distance f* cdt'/a(t') and the other quantities in ( JI3"]) as having assumed their 



3 This is a consequence of the Wiener-Khinchine theorem and the fact that 
a dispersive medium does not change the energy spectrum of the radiation, see 
e.g. http : //light . ece . illinois . edu/ECE460/PDF/LCI .pdf . 
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values at time t. This assumes of course that the evolution of n e is caused solely by the 
expansion. For sources out to at least z — 1, the warm IGM indeed maintains its constant 
baryonic mass fraction of ~ 50 %, i.e. additional (evolutio nary) corrections are not expected 



to be important; see Figure 2 of ICen and Ostriker I (119991 ) 



Two features are apparent in (IT31 . First, in calculating the default value of ( IT3|) we 
assumed Auo/ujq = 0.01, which in the context of astrophysical sources is the spectral width 
of a narrow emission line or the telescope frequency selection. Second, we list in Table 1 
the column density of the various plasma components that extragalactic light passes through 
before reaching a ground based telescope; it is evident that apart from directions intercepting 
rich clusters of galaxies the IGM delivers the largest effect, which is why we expressed ([TBI 
in representative units of n e for this component. 

Turning to observational strategy, we shall focus upon distant quasars, i. e. point sources, 
to ensure spatial coherence in the light. From Fig. 3 of the 5,000 quasar sample of 



Singal et al I ( l2012h one may assume the radio spectral index of 0.6 to obtain the intrinsic 
luminosity density of 10 31 ergs s _1 Hz -1 at 1 GHz as a representative estimate for quasars 
with z < 0.25 (£ < 1 Gpc). Applying a gaussian spectral filter of Au/uq = 0.01 to the 
arriving continuum radiation (the spectrum of which is taken as flat), one obtains r c = 1.6 
xl0~ 8 s, f « 5 x 10 5 from ([13]), and a = 8 x 10 3 s" 2 and b = 4 x 10 9 s~ 2 from ©. The 
photon 'occupation number' will be n « 20 if a quasar at 1 Gpc distance is observed by a 
telescope of diameter 300 m (Arecibo, it should also be mentioned that over the beam of this 
telescope the equivalent number of cosmic background photons is n CMB « 3), and indicates 
that one is in the classical (phase) noise limit for the radiation intensity, see the previous 
section. An exposure time of 1 hour would yield N rj 4 x 10 12 as the number of collected 
photons. 

Thus the numbers all point to the validity of the classical description of radiation given 
in the previous section, which also presented a relatively straightforward means of using the 
quasar light to measure line-of-sight dispersion, viz. the microscopic random fluctuations 
should have their highest amplitude on the timescaleS a -1 / 2 ~ 10 millisecond and not the 
(much shorter) coherence time r c , although the Fourier transform of the light curve should 
reveal fluctuation power on all scales between a -1 / 2 and r c , monotonically falling from the 
former to the latter. By observationally determining a _1//2 in this way, the line-of-sight IGM 
plasma column may then be inferred from (E^i), ( fl3|) . and knowledge of the bandwidth Au 



4 In addition to having a sufficient number of photons per such an interval, radio observations at this timing 
resolution is also not expected to be a problem, since the highest resolution ever reached is 4 nanosecond 



(jHankins et al I (|2003f )) 
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of frequency selection. The crucial datum here is the upper cutoff timescale of a -1 / 2 , i.e. one 
does not expect any significant variations in the intensity above this scale, until the scale ~ 1 
day when genuine changes in the emission properties of the quasar source have been noted by 
e.g. lEggers et al I (120001 ). To the best of the authors' awareness no attempt have been made 



to monitor quasars on the millisecond time frame. Of course, if source variability occurs 
down at this range the proposed technique will fail. It is extremely difficult to envisage such 
phenomena, however, because of the size of the quasar emission region. 
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Type of plasma 


n e (in cm 3 ) 


Column length I 


Column density n e £ (cm 2 ) 


IGM 


10^/10.7 


3 Gpc 


9 x 10 20 


Rich clusters of galaxies 


10~ 3 


2 Mpc 


6 x 10 21 


Interstellar medium (ISM) 


0.03 


1 kpc 


9 x 10 19 


Interplanetary medium (IPM) 


10 


100 AU 


1.5 x 10 16 


Earth's ionosphere 


10 5 


300 km 


3 x 10 12 



Table 1: Properties of plasmas on various length scales, with the following noteworthy points, 
On the ISM (interstellar medium) the reader should beware of the anisotropy in n e presented 
by the Galactic disk, i.e. the column density as t abulated remains representative of any 
direction including high Galactic latitudes (see e.g. iHowk et al I (120061 ) ) . assuming that one 
avoids the disk when performing cosmological observations. On the IGM, the density n e is 
obtained by assuming that the baryonic IGM consists principally of the 10 5-7 K plasma (the 
WHIM, or warm-hot intergalactic medium) of sectio n 1, with normalized cos mic density 
^whim ~ 0.02 between z = and at least z — 1 (ICen and Ostriker I (119991 ) . Fig. 2b), 
where h , 7 denotes the present Hubble constant H in units of 70 km s" 1 Mpc -1 and £ is the 
comoving distance J dt/a(t). 



